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Abstract 

We calculate the scalar semileptonic kaon decay in finite volume at the momentum 
transfer tm = {rriK — "i-tt)^, using chiral perturbation theory. At first we obtain the 
hadronic matrix element to be calculated in finite volume. We then evaluate the finite 
size effects for two volumes with L = 1.83 fm and L = 2.73 fm and find that the 
difference between the finite volume corrections of the two volumes are larger than 
the difference as quoted in [18] . It appears then that the pion masses used for the 
scalar form factor in ChPT are large which result in large finite volume corrections. 
If appropriate values for pion mass are used, we believe that the finite size effects 
estimated in this paper can be useful for Lattice data to extrapolate at large lattice 
size. 



PACS: 11.15.Ha, 12.39.Fe, 13.20.Eb, 14.40.Aq. 



1 Introduction 



A precise determination of the CKM matrix element, Vus, is highly demanded if one expects 
to find a footprint of new physics in the unitarity relation of the first row of the CKM matrix 
[H [2] , namely, 

|Kd|' + |KsP + |K6|' = i. (1) 

Among three entities above, \Vus\ is less precisely known from experimental measurements 
in strangeness changing weak kaon decays K13 and therefore dominates the uncertainty in 
the unitarity relation. The theoretical framework for the determination of this quantity 
is explained by Leutwyler and Roos in [3j. To achieve a more precise value of |V^<i|, one 
requires a better control over the long distance effects due to the strong interactions in the 
process under consideration. In fact one can only obtain experimentally, the combination 
|\4s/+(0)| by measuring the semileptonic kaon decay rates. Here, /+(0) denotes the related 
vector form factor at zero momentum transfer. This can be exemplified in the master 
formula for K13, see [1] and references therein 

r {k^ n^t.,) = C-'^^^^ \v^sfr'''\n' (1 + 2A^i,,) , (2) 

where Gp is the fermi constant. Sew incorporates the short distance electroweak correc- 
tions and I^/ is the phase-space integral depending on the slope and curvature of the form 
factors. /+(0) contains the non-perturbative effects of strong interactions whose precision 
estimate may lower the uncertainties on \Vus\- There exist two main approaches for the 
calculation of the vector form factors. One powerful technique is the application of chiral 
perturbation theory (ChPT) as an effective field theory describing low energy strong in- 
teractions. The determination of /+(0) at one loop order within ChPT was evaluated by 
Gasser and Leutwyler [5]. Along the same lines the calculations were extended further to 
incorporate two loop order corrections Despite its model-independent prediction, how- 
ever, ChPT calculations may become less predictive at higher orders due to undetermined 
parameters. 

Lattice QCD on the other hand, is a brute force approach to determine the QCD 
observables by evaluating the functional integral of QCD numerically. In recent years the 
technical procedures have reached to a point that calculations with relatively small quark 
masses are feasible, see [3 El IHl UHl IH] and references therein. Lattice data are not free 
from systematic errors e.g. finite volume effects and chiral extrapolation. As soon as 
pion mass is small enough and lattice size is large enough in order to make legitimate 
application of ChPT for the extrapolations, one can get reliable predictions out of Lattice 
QCD data. There are a large amount of works in the literature on exploiting ChPT to 
correct the systematic errors. Lattice size effects on pion mass are studied in [T2] and for 
pion decay constant in [131 [II]- Along the same arguments this type of calculations are 
done for quark vacuum expectation values [15]. Moreover, finite volume effects on pion 
pion scattering parameters are evaluated near threshold [16]. There exist Lattice data 
on the semileptonic vector and scalar form factors, see for example [171 [TSl HSl [201 1^ - 
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As suggested in [TTl [22] the starting point in doing these calculations in the Lattice is 
to obtain the scalar form factor at the maximum momentum transfer tm = {itik ~ ""^tt)^, 
where high statistical accuracy is achievable. We then deem interesting to evaluate the 
finite size effects for the scalar semileptonic kaon decay. This decay gives the important 
key role in kaon physics. Therefore we set the main motivation of the present work as 
pushing the line of the ChPT applications to a case which involves external momentum. 
For the sake of simplicity we assume external momentum takes on zero spatial component. 
To this end, we have checked how far we can get on extending the reliability of ChPT to 
doing chiral and finite volume extrapolations. Earlier work on finite volume calculations 
for meson current matrix element with external spatial momentum can be found in [23] . 

The outline of this article is as follows. We begin with an introduction about ChPT in 
finite volume and infinite volume in Sec.|2j On Sec. [3] the definition of the Ki^ form factors 
are reviewed and in Sec. H] one loop result of the hadronic matrix element in infinite volume 
and in a tensor form is provided. Finite volume calculations entailing necessary Feynman 
integrals in finite volume are sketched in SecEJ The last section reports on the finite size 
effects on the scalar form factor and ends with conclusion. 

2 ChPT in infinite and finite volume 
2.1 ChPT in infinite volume 

Chiral perturbation theory (ChPT) is an effective field theory which describes the low 
energy dynamics of Quantum Chromo Dynamics (QCD). Weinberg in a seminal paper 
for the first time paved the way for a systematic application of effective Lagrangian |24] . 
Pseudo-Goldston mesons play the role of dynamical degrees of freedom in the effective 
Lagrangian. These mesons are resulted as a consequence of the spontaneous breakdown 
of the chiral symmetry in QCD. The expansion parameter is generically in terms of the 
momentum, p^, and quark mass, ntq. Later on Gasser and Leutwyler in two elegant papers 
extended the Lagrangian to order |25l |26] . One can arrange the full effective Lagrangian 
as series of operators with increasing number of derivatives and increasing power of quark 
mass 

Cff = A + A + A + ■ ■ ■ • (3) 

Here the subscripts show the chiral order. The lowest order Lagrangian involving two terms 
has the following form 

A = ^ {{D.UD^U^) + (xf/^ + Ux^)) . (4) 

The mass relation = Bq {niu + rrid) allows one to count quark masses as order p"^. The 
next to leading order Lagrangian consists of 10+2 independent operators with correspond- 
ing effective low energy constants (LEC's) 

£4 = Li{D^U^D^'Uf + L2{D^U^D,U){D^U^D'U) 
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+L3{D''U''D^UD''U^DM) 
+L,{D'^U^D^U){x^U + xU^) 
+L,{D''U''D^U{x^U + C/tx)) + Leix^'U + xU^f 
+Lr{x''U - xU^f + Lsix^'Ux^'U + xU^U"^) 

+H2{X^X) ■ 



(5) 



The notation (...) = Tip (...) indicates the trace over the flavors. The matrix U e SU{3) 
parameterizes the eight hght mesons in an exponential representation 



U{(j)) = exp(iV2(/)/Fo) 



(6) 



where 



V2 ye 



TT ' 



V 



TT 



K- 



7^3 ^ ^8 

V2 \/6 

^0 



(7) 



The covariant derivative and the field strength tensor are defined as 

D„U = dJJ - irJJ + iUl^ , Fl:, = d^l, - dj^ - i [l^, I,] 



(8) 



and a similar definition apply to the right-handed field strength. Here and represent 
the left-handed and right-handed chiral currents respectively, x is parameterized in terms 
of scalar (s) and pseudo scalar {p) external densities as x = 2Bq {s -\- ip). For the process 
discussed in this article it is sufficient to set 



m 



m 



= 0. 



92 

V2 



V v*w- 

\ us a 



(9) 



The weak coupling constant, 5f2, is related to the Fermi constant by 

8m^ ~ y/2' 



(10) 
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2.2 ChPT in finite volume 



The effective Lagrangian approach is also apphcable to a system which is bound to a large 
space of volume V = with an infinite time extent. The effective Lagrangian to be used 
for a system enclosed into a finite size is the same as one we use for processes which occur in 
infinite volume with the same values for the low energy effective constants. The pioneering 
works [271 12HI [29] provide the foundations in this regards. To do the calculations in finite 
volume as suggested in |29] one imposes a periodic boundary condition on particle fields 
which results in the momentum quantization and consequently the modifications of the 
quantum corrections. The propagation of a particle in space-time is defined through the 
two-point correlation function. Therefore the correlation function in finite volume becomes 



Gv = ^,T.J^Giq'^q), (11) 



where the three dimensional momentum takes the discrete values 

^ 27r^ 

q= —n, (12) 

with n is a three dimensional vector with integral components and L is the linear size of 
the box. Since the reliability of ChPT is limited to the small momenta, in finite volume 
this leads to the condition that 

F^L»1, (13) 

where is the pion decay constant. In addition we stay in a limit where the Compton 
wavelength of the pion is smaller than the lattice size, L, so that the zero mode of the pion 
field does not become strongly coupled. This condition is satisfied if 

m^L » 1 . (14) 

This is the so-called p-regime in which we do our calculations here. Moreover, the quantity 
m^L plays the role of the power counting in the perturbative calculations. 

3 The definition of the Kis form factors 

The semileptonic kaon weak decays which are generally shown as Ki^ are: 

^+(p)^vr°(p')/+(p/)^/(p.), (15) 
K\p)^iv-{p')t{pi)ui{p,). (16) 

Here / stands for electron and muon. The two other processes are charge conjugate modes 
of the processes above. 



Neglecting scalar and tensor contributions, the matrix element takes on the structure 

m = ^v:,j^f;:{p\p), 
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with 



F^{p',p) = < 7c%p')\s^,umK+ip) > . (18) 
The general form of the vector matrix element is defined 

<7r%p')\rs^,u{0)\K^{p) = 



The matrix element for neutral kaon can be obtained when we replace by 



(19) 



F-{p',p) = < n-{p')\s^,umK\p) > 



(20) 

The four Ki^ form factors f^~^'"°{t) and f^°'^ (t) depend on the four momentum squared 
which is transferred to the leptons 

t = {p'-pf = {pi+Puf. (21) 

In the isospin limit f^^'^°(t) = f^°'^ (t). In this article the matrix element in infinite 
volume is obtained for the neutral koan decay. The scalar form factor which describes the 
S-wave projection of the matrix element is defined as 

m = Mt)+ , ^ j .it). (22) 

At the point of zero recoil, the general decomposition of the vector current matrix element 
becomes 



<T,-{p')\s^,umK\p) >t. 

^ ■.[{rriK + m^)ff^'{t^) + (m^ - m^)f'^"^'{t^)]6^o , 



V2' 

(23) 

where tm = {ttik — "n^Tr)^ is the value of momentum transfer at the special point in the 
kinematical region. Given the definition of the scalar form factor in Eq. (!22|) and the matrix 
element in Eq. ( 123|) we find the following relation at the maximum momentum transfer 



(24) 
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(a) (b) (c) (d) (e) 

Figure 1: Feynman diagrams contributing to the form factors. Diagram (a) is of 0{p'^) with 
a vertex from the 0{p^) Lagrangian. The rest of the diagrams are of 0{p'^). Diagram (d) 
takes its vertex from the 0{p^) Lagrangian and diagram (e) is the wave function corrections 
to the external legs. The thick lines define meson propagations and wavy lines indicate the 
gauge boson. 



4 Analytical results in infinite volume 

In this section we recalculate the Ki^ vector form factors at order 0{p'^) and 0{p'^) in 
the isospin limit. On top of this, we also obtain a new expression at one loop order for 
strangeness changing matrix element in a tensor form which is useful for finite volume 
calculations to be discussed in Sec. [5] and Sec. [61 In the next subsection we introduce the 
matrix element in its new form and recall the old results. Feynman diagrams up to order 
in chiral order needed for the interested process are provided by Fig. [H 



4.1 0{p^) and 

In the isospin limit the vector form factors and scalar form factor at leading order are 

fl{t) = m = l, f_{t) = 0. (25) 

The superscribe indicates the chiral order. These results are obtained by calculating the 
tree Feynman diagram in Fig. [T]^a). At the next order there are four diagrams contributing 
to the matrix element as depicted in Fig. [T](b-e). Without using the tensor simplification 
relations for the tensor integrals, we evaluate all the diagrams and sum of the four diagrams 
read the following expression 

F{p,p').e = j-,[2q'L, + [^A{ml) + ^-A{ml) 

+ ^B,,{ml, m% q^)Ye^ + [-2(m^ - ml)U 
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1 

2' 



2\ rD/2 2 2 



12 

2 _2 2 



^6 



-^g'e'^] + B^{m\, m% q^)^^{p + p'fq.e 
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Abb 



-qh^] 



(26) 



with 



p' + p, 



q = p -p, 



(27) 



and e is the polarization four vector of the gauge boson. The integral functions A,B,Bf^ 
and B^i, are introduced in the Appendix. In infinite volume, it is possible to use the 
tensor simplification relations which are provided in the Appendix, to reduce the tensor 
integrals in terms of scalar functions. The one loop integrals can be evaluated by applying 
the dimensional regularization scheme. The infinities arising from loop integrals can be 
canceled by redefining (in a renormalization scheme as discussed in [26]) the LEC's in 
terms of renormalized LEC's and subtracted parts containing infinities. After doing the 
renormalization, we ensure the cancellation of infinities and recover the old results [3] 



3^, 



3^. 



TT 



2L;t + ^A{ml) + -A{ml) 



2. 3- 
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--B22{ml,ml,t) 



(28) 



f. 
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+2mlLl-^A{ml) + l^A{ml) 
+\ A{ml) + B{ml, m^, t){^t- 
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^2 "^^) + B{m\, ml, t){^t-^m\ 
+-^m^)+ Bi{ml,m\,t){-^t + ^m\ 
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5 Finite volume calculations 

When we enclose a system in a finite box, Lorentz symmetry is no longer respected. This 
can make it impossible to define without ambiguity vector form factors for such a system. 
It also turns out impossible to define vector form factors at the point of zero recoil energy 
in finite volume because it is always possible to reshuffle terms in the matrix element 
in such a way to obtain different values for /+ and /_ according to their definitions in 
Eg. (1231) . However, we see from Eq. llMl) that the time component of the matrix element 
at the maximum momentum transfer, tm, is proportional to the scalar form factor and 
therefore the evaluation of the time component of the matrix element in finite volume can 
give us the scalar form factor at t^. The one- loop Feynman integrals which contribute to 
the matrix element of the process are evaluated in this section. We expect finite volume 
effects arising from the simultaneous propagation of kaon and eta particles to be much 
smaller than the same effects when pion and kaon propagate simultaneously in the loop. 
We demonstrate this claim quantitatively in the next subsection. 



5.1 Scalar one-loop integrals 

We work out in this subsection the one-loop scalar integrals in detail. In what follows we 
define for a generic function, F, AF = Fy — F^o- Subscripts V and oo indicate integra- 
tion in finite and infinite volume, respectively. In finite volume, as we mentioned before, 
momentum q gets discrete values, q = ^n, but we obtain our formula for the case with 
infinite time extension. We begin with the tadpole integral 



dqo ^ f d^q 



To get the second equality we have used the Poisson summation formula 

^ E f{^') = I + E / i^sfi^iy''-'- (31) 



iLq.n 



1= 



8 



By evaluating the contour integral over go in Eg. (1301) we arrive at 



^iLq.n 



We finally take the final integral and obtain 

AA = Y: j^Mk)K,{mL\n\), (33) 

where function Ki is the modified Bessel function of order one. m(k) accounts for the 
number of ways that the equality = + n| + rzg is satisfied for positive and negative 
integer numbers for rii, n2 and 71,3. This result agrees with that in [T5] . 

The second Feynman integral we need to consider here, involves the propagation of two 
particles with different masses. At the momentum transfer, tm = [Mk — M^r)^, we first 
perform the integral for particles with masses M = Mk and m = in the propagators 



(34) 

where the Poisson summation formula is employed to get the second line. For a simplifying 
assumption we take the external momentum completely in the time direction as p = {po, 0) 
and we proceed by exploiting the Feynman parameter formula to get 

1 rdqo ^ f d^q e^^^'" 



Bv = Boo - i dx —Y^ 



27r (27r)3 [g2 + x{l - x)pl - xm? - (1 - x)M' 



212 



(35) 



We first perform a contour integration over go and then plug in = to find 

4^o7 (27r)3 Jo [g^ + (a;m + (a;-l)M)2]3/2- ^^'^^ 
Taking the integral over x gives 

A{M -m) {2rr)^ ^cP^q' + M'^ q^^^Tm^ 

The second piece contains a three dimensional integral where integration over the angular 
part leaves us with 

i ^ 1 /■+00 Me^\^\^\ me'^\^\^\ . 

^v = ^oo- TT^TTTT, E ^ / ^1^1 



27r)24(M-m)L 7-00 \q\^\q\^ + \q\\l \q\^ + m'' 

{31 
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To carry out the last integration we make use of the convolution technique. After taking 
the integral we obtain 



= TTTTiT ^ E M[L^,{ML\n\)Ko{ML\n\) + Lo{ML\n\)K,{ML\n\)] 

—m[L_i{mL\n\)KQ{mL\n\) + LQ{mL\n\) Ki{mL\n\)]. 

(39) 

La and are respectively modified Struve and modified Bessel functions of order a. In 
Fig. |2l the function AB is plotted versus L for three different values of the pion mass while 
we change the kaon mass correspondingly. The results in Fig. |2] show a reasonable trend. 
In addition, we need to perform another one-loop integral for propagating particles with 
masses and Mk- It turned out that one cannot find a formula in closed form in this 
case. The final result for the K — rj loops reads 

^Q^J_Sr ^ (H1 + H2 Sin{L\n\q)qdq 



87r2^^oL|n|V Jo + H3)^{Hi + H^f + H3 + 

Sin{L\n\q)qdq 



(40) 



where 



H^ = Mk-M^, 

_ Ml - - {Mk - M^f 
' ~ 2{Mk - M^) 

H, = Ml - Hi (41) 

We solve the integral above numerically and plot the function Ai? versus L as depicted in 
Fig. 131 By comparing the results in Fig. |2] and Fig. |3l we conclude that the K — rj loops 
are suppressed with respect to the n — K loops. 



5.2 Tensor one-loop integrals 

The tensor integral defined in the Appendix has only one component, i.e. -B°, needed 
to be evaluated in the finite volume at the momentum transfer tm- We thus consider the 
integral 

_ i rdqo^ go rO , /" /" d^q 



L^J 27r ^ {q^-m^){{q + py-M^) °° 'J 2ti i^J (27r 

„ ALq.n 

X (A2) 

{q^ ~m?){{q + pY - M-^y ^ ' 
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Figure 2: tt-K loops in finite volume for three different values of the pion mass and kaon 
mass. 



The second term which we call AB'^ can be written in terms of the functions which we 
obtained already. After doing some manipulations on the integrand, we obtain 

Afi° = — — i -[AA(m) - AA(M) + 2m(M - m)AB]. (43) 

2(M — m) 

The last integral we need to calculate in the finite volume is the B^'^ component of 
the tensor integral B'^'^. Again we apply the Poisson summation formula to arrive at the 
second line for the integral below 

ooo =_A/"^V ^ 

^ J 27^ ^ {q" - m?){{q + pY - M^) 

= B00_^ fdqoy ^° (44) 

7 27r £^^7 (27r)3(g2-m2)((g + p)2-M2)' ^ ^ 

We integrate over go following the same argument sketched before and find 



(27r)3 2Vg' + M'' A{M-m)J (27r)3 
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Kaon-Eta Loops 




Figure 3: K-rj loops in finite volume for three different values of the pion mass and kaon 
mass 



Substituting the results of Eg. ( 15^ and Eq.( lS7|l into the relation above will yield us 

AS°° = AA{M) - — — -\MAA{M) - mAA{m)] + m'^AB. (46) 

2(M — m) 

5.3 An asymptotic formula for the scalar form factor 

It is important to know the behavior of the finite volume correction of the scalar form 
factor in asymptotically large volume. To this end, we provide a formula in which we have 
ignored eta loop effects as well as terms with kaon mass in the exponential and obtain the 
following result 



AfoiU, 



1 / y/2m^{mK — 19m.n-) 



F2V24(27rL)3/2(m, + mi^) 
E -TZ^- (47) 



\n\ 



Hence, one clearly notices that the volume corrections to the scalar form factor are expo- 
nentially suppressed. 
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1.06 - 



1.04 



1.02 



m'k = 500 MeV 
mk = 600 MeV 




0.45 



Figure 4: We compare the pion mass dependence of the scalar form factor at the endpoint 
in the infinite volume for two different values of the kaon mass. 



6 Numerical results and conclusions 



It is seen that the value of the scalar form factor at the endpoint in infinite volume only 
depends on one low energy constant i.e., LI, which we use as input = 0.97205 x 10~^, 
see |30] for detailed discussion on this LEG. The subtraction scale is /i = 770 MeV. There 
is an ambiguity on the value we should put in for the pion decay constant when one 
calculates a quantity at one loop order. We use the physical value for the pion decay 
constant instead of its lowest order value since the difference shifts to higher order. The 
value = 0.092AGeV is chosen as input. The loop integrals in infinite volume are 
performed with the renormalization scale n = 770 MeV. We show in Fig. Hlthe pion mass 
dependence of the scalar form factor in infinite volume at the endpoint for rriK = 500 MeV 
and rrix = 600 MeV. Moreover, from Fig. HJ one can see that in going away from the 
physical value of the pion mass to an unphysical value of about 350 MeV the scalar form 
factor in the infinite volume changes about 5 percent. In addition, we study the volume 
dependence of the form factor at order by defining the ratio 



R 



fo {tm) /o°(^m) 
fo°(tm) 



(4J 



The ratio Rj^ is plotted in Fig. [5] versus lattice size, L, for three different values of the 

300 MeV and = 450 MeV, with a fixed value 



pion mass, i.e. m,, 



220 MeV, 
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Figure 5: The ratio Rf^^ is plotted versus the hnear size of the lattice for pion masses 
= 220 MeV, m-,^ = 300 MeV and = 450 MeV for a fixed value of the kaon mass, 
niK = 500 MeV. 
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m^(GeV) 


niKiGeV) 




A/o (Only n — K loops) 


A/o (Including K — r) loops) 


Jo 


/,r :t8: 


0.329 


0.575 


4.50 


0.08459 


0.08490 


1.09532 


1.02143(132) 


0.416 


0.604 


5.69 


0.11169 


0.11196 


1.11360 


1.00887(89) 


0.556 


0.663 


7.61 


0.12323 


0.12340 


1.12179 


1.00192(34) 


0.671 


0.719 


9.19 


0.11435 


0.11446 


1.11382 


1.00029(6) 



Table 1: We compare our finding for tlie scalar form factor in finite volume with those 
quoted in [18], corresponding to the pion masses and kaon masses given in the first and 
second column. The linear size of the lattice in these calculations is 2.74 fm. Finite volume 
corrections to the scalar form factor without and with K-rj loop contributions are given in 
the forth and fifth column respectively. 



for the kaon mass, mx = 500 MeV. The results in Fig. [5] clearly show that for lattice 
sizes smaller than L = 6.15 fm there is a larger finite size effects for the larger pion mass 
while the ratio behaves differently for the lattice sizes larger than L = 6.15 fm. We found 
out that the term with the function in our finite volume expression for the scalar form 
factor, takes responsibility for finite volume effects to become large for larger pion mass at 
small volumes. The reason is that the function B'^ in the finite volume does not converge 
fast enough to compensate the quadratic growth of the pion mass in front of B^. However, 
as we see from Eq. HTl the finite volume corrections to the scalar form factor exhibit a 
standard behavior at asymptotically large pion mass, wherein the corrections suppress 
exponentially. We therefore stress that the finite volume effects for the scalar form factor 
at the endpoint can be sizable for L = 4 fm for large pion mass of about 300 MeV. 

In Table 1. we present the finite volume corrections to the scalar form factor and the 
value of the form factor in finite volume with L = 2.73 fm for particular values of the 
pion mass and kaon mass as used in [TS]. Given the systematic errors in the scalar form 
factor quoted in [18] our results in Table. [T]show that the lattice data is precise enough for 
heavier pion masses to take into account the r] contributions at finite volume. However, our 
results for the finite volume correction to the scalar form factor at finite volume appear to 
be large so that the resulting scalar form factor in finite volume is larger than those quoted 
in [18] while for light pion masses the discrepancy becomes smaller. This is expected from 
Fig. O since our calculated finite volume corrections become large for large pion masses 
while being in the p-regime, and therefore this is the question of how much SU(3) ChPT is 
reliable for larger pion masses in the process studied in this article. Moreover, we have also 
evaluated the finite size effects for a smaller lattice size with L = 1.83 fm and find that 
the finite volume effects here is about 40 percent which exceed the almost 10 percent finite 
volume corrections obtained for a lattice with L = 2.73 fm. The difference between the 
estimated finite size effects from ChPT for L = 2.73 fm and L = 1.83 fm is quite larger 
than the amount of difference as quoted in Table III found in [18]. This is again due to the 
fact the used pion masses in our ChPT calculations for the scalar form factor are large. 

Further work in this direction is to extend our results to non-vanishing spatial momen- 
tum since lattice practitioners use these values to interpolate to zero momentum transfer. 
Moreover, our study can be extended to the case of twisted boundary condition. 
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7 Appendix 

We introduce the one loop Feynman integrals as follows 

A^m") - - / ,5^, / ^ , (49) 



B,,{m\M\p') =-[ 

% J 



{q^ -w?){{q + pf - W^y 
1 r df^q 



(51) 



d 



i J (27r) 

{q''-'m?y{q + py-M-'y ^ ' 

Invoking lorentz symmetry, the tensor integrals above can be written in terms of scalar 
functions 

B^{m\M\p')=p^B^{m\M\p'), (53) 

B^,{m\M\p') =p^p,B,2{m\M\p') 

+g^.,B22{m\M\p'). (54) 
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